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ABSTRACT 

Reverberation mapping of nearby active galactic nuclei has led to estimates of broad-line-region (BLR) sizes 
and central-object masses for some 37 objects to date. However, successful reverberation mapping has yet to be 
performed for quasars of either high luminosity (above L opt ~ 10 46 erg s -1 ) or high redshift (z > 0.3). Over the 
past six years, we have carried out, at the Hobby-Eberly Telescope, rest-frame-ultraviolet spectrophotometric 



monitoring of a sample of six quasars at redshifts z = 2.2-3.2, with luminosities of L opt ~ 10 404 — 10 erg s~ 
an order of magnitude greater than those of previously mapped quasars. The six quasars, together with an 
additional five having similar redshift and luminosity properties, were monitored photometrically at the Wise 
Observatory during the past decade. All 1 1 quasars monitored show significant continuum variations of order 
10-70%. This is about a factor of two smaller variability than for lower luminosity quasars monitored over 
the same rest-frame period. In the six objects which have been spectrophotometrically monitored, significant 
variability is detected in the C IVA1550 broad emission line. In several cases the variations track the continuum 
variations in the same quasar, with amplitudes comparable to, or even greater than, those of the corresponding 
continua. In contrast, no significant Lya variability is detected in any of the four objects in which it was 
observed. Thus, UV lines may have different variability trends in high-luminosity and low-luminosity AGNs. 
For one quasar, S5 0836+71 at z = 2.172, we measure a tentative delay of 595 days between C IV and UV- 
continuum variations, corresponding to a rest-frame delay of 188 days and a central black-hole mass of 
2.6xlO 9 M . 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert — Quasars: general 
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1. INTRODUCTION 

Reverberation mapping has been used in the past two 
decades to estimate the sizes of the broad-line-emitting re- 
gions (BLRs) in several dozens of active galactic nuclei 
(AGNs), and thus to infer the masses of the black holes at their 
centers (e.g., Kaspi et al. 2000). The technique utilizes the 
light-travel-time delayed flux response of the BLR to changes 
in the continuum flux (for reviews, see Peterson 1993; Netzer 
& Peterson 1997; Peterson 2006 and references therein). Re- 
cently, a compilation of all available reverberation-mapping 
data, for 37 AGNs, was analyzed in a uniform and self- 
consistent manner. The relationship between luminosity, L, 
and the BLR size, Rblr, in several luminosity bands was stud- 
ied by Kaspi et al. (2005), and the black-hole mass-luminosity 
relation, Mqh-L, based on these data was re-derived by Pe- 
terson et al. (2004). Current reverberation studies cover four 
orders of magnitude in AGN luminosity in a well-sampled 
manner, from ALa(5100) ~ 10 42 ergs s -1 to 10 46 ergs s -1 . 
In this range, R^lr oc L a , with a = 0.67 ±0.05 for the op- 
tical continuum, confirming the relations found by Kaspi et 
al. (2000; c.f., Bentz et al. 2006 find a = 0.518 ±0.039) and 
M B h oc L , with (3 = 0.79 ± 0.09. In a recent study Peterson 
et al. (2005) measured the BLR size of the C IVA1550 emis- 
sion line in the least-luminous AGN, NGC 4395, which has an 
optical luminosity of ALa(5100 A)=5.9x 10 39 ergs s -1 . They 
find the C IV BLR size to be 1 ± 0.3 light hr, consistent with 
the size expected from extrapolating the Rblr-L relation to 
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lower luminosities. 

The relations between AGN luminosity, BLR size, and 
black-hole mass have been widely used for "single-epoch" es- 
timates of black-hole masses and accretion rates in studies ad- 
dressing the issues of AGN accretion history and black-hole 
growth (e.g., McLure & Dunlop 2004; Barger et al. 2005; Yu 
et al. 2005; Kollmeier et al. 2006). However, these studies 
necessarily extrapolate the relations to luminosities and red- 
shifts well beyond the ranges in which they were measured 
(e.g., McLure & Jarvis 2002; Vestergaard 2002; Woo & Urry 
2002; Vestergaard & Peterson 2006), since the only measured 
relations are in the luminosity range 10 42 -10 46 ergs s" 1 and 
for redshift < 0.3. As a result, the conclusions of such studies 
depend on the untested assumption that these extrapolations 
are valid. Although a posteriori explanations of the physi- 
cal plausibility of the observed relations can be found, it is 
quite possible that subtle or strong deviations from the rela- 
tions occur at high luminosities or redshifts (Netzer 2003). It 
is well known that some other AGN properties, such as the 
optical-to-X-ray spectral slope and emission-line equivalent 
widths, depend on luminosity (e.g., Baldwin 1977; Strateva et 
al. 2005; Steffen et al. 2006). AGN outflows also seem to be 
different between Seyfert galaxies (several hundred kms" 1 ) 
and broad absorption line (BAL) quasars (several thousand to 
tens of thousands km s -1 ), and thus might depend on luminos- 
ity (Laor & Brandt 2002). Actual reverberation measurements 
for AGN with high luminosities and redshifts are therefore de- 
sirable. 

However, reverberation mapping of high-luminosity 
quasars is an ambitious undertaking. Quasars of the highest 
luminosities (with bolometric luminosity, Lbol ~ 10 47 — 10 48 
ergs s ) are expected to harbor some of the most massive 
black holes known, with Mbh ^ 10 9 M Q . More massive 
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TABLE 1 
Object Characteristics 



Object 


RA (2000) 


Dec (2000) 


my 


Redshift 




A/spec 


AL A (5100A) 




-^^comp 


R 


Exp. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


CO 


(8) 


(9) 


(10) 


(11) 


(12) 


S4 0636+68 


6:42:04.2 


67:58:36 


16.6 


3.180 


90 


11 


47.28 


131.9 


36.6 


133.2 


900 


S5 0836+71 


8:41:24.3 


70:53:42 


16.5 


2.172 


70 


16 


46.81 


66.7 


341.2 


10064.5 


600 


SBS 1116+603 


11:19:14.3 


60:04:57 


17.5 


2.646 


85 


15 


46.92 


147.2 


62.7 


632.7 


900 


SBS 1233+594 


12:35:49.5 


59:10:27 


16.5 


2.824 


76 


15 


46.97 


125.6 


79.0 


1.3 


600 


SBS 1425+606 


14:26:56.2 


60:25:51 


16.5 


3.192 


90 


21 


47.43 


78.6 


2.4 


4.6 


300 


HS 1700+6416 


17:01:00.6 


64:12:09 


16.1 


2.736 


96 


17 


47.49 


91.4 


135.4 


3.6 


300 


S5 0014+81 


0:17:08.5 


81:35:08 


16.5 


3.366 


82 




47.56 






493.3 




S5 0153+74 


1:57:34.9 


74:42:43 


16.0 


2.338 


65 




46.85 






12377.8 




TB 0933+733 


9:37:48.8 


73:01:58 


17.3 


2.528 


86 




47.07 






5.2 




HS 1946+7658 


19:44:55.0 


77:05:52 


15.8 


2.994 


102 




47.63 






1.7 




S5 2017+74 


20:17:13.1 


74:40:48 


18.1 


2.187 


77 




46.44 






2765.8 





NOTE. — Upper entries are objects monitored both photometrically and spectrophotometricaily. Lower entries are objects monitored only photomet- 
rically. Col. (1) Object name. Cols. (2)-(3): Right ascension and declination (2000) from NED. Col. (4): V magnitude from the Veron-Cetty & Veron 
(1993) Catalogue. Col. (5): Redshift from NED. Col. (6)-(7): Number of photometric and spectrophotometric observing epochs. Col. (8): Log luminosity 
at 5 100 A (luminosity in units of ergs s~') computed from the rest-frame flux density at 1450 A of the averaged spectrum, using a power law of f v oc f~°' 5 , 
and correcting for Galactic absorption. Cols. (9)-(10): Angular separation in arcseconds and position angle in degrees of comparison star (see § 12. 31 . Col. 
(1 1): Radio (1.4 GHz)-to-optical (estimated at 4400 A) flux ratio. Col. (12): HET/LRS exposure time in seconds (two consecutive exposures with this 
exposure time were obtained each visit). 



black holes may have slower continuum flux variations with 
smaller amplitudes. This has been observed in X-rays (e.g., 
Lawrence & Papadakis 1993; Uttley et al. 2002; Markowitz et 
al. 2003; O'Neill et al. 2005), and possibly in the UV-optical 
bands (e.g., Giveon et al. 1999; Cid Fernandes et al. 2000; 
Vanden Berk et al. 2004). The required observing periods of 
high-luminosity quasars are also significantly lengthened by 
cosmological time dilation, since such sources are typically 
found at high redshifts (z > 1). On the other hand, the ability 
to monitor high-z objects in the rest-frame UV, in which AGN 
variability amplitudes are routinely higher than in the optical, 
can lead to better characterized continuum light curves. The 
smaller intrinsic variability amplitude of the continuum could 
result in smaller flux-variability amplitudes for the emission 
lines, affecting the ability to detect the time delay in the 
BLR response. Furthermore, since the equivalent widths 
of high-ionization emission lines decrease with increasing 
luminosity (the well-known 'Baldwin Effect'; Baldwin 
1977), fluxes of emission lines such as C IVA1550, which 
are most suitable for such studies, are harder to measure. 
On the other hand, the fluxes in the UV lines are generally 
higher than those of the Balmer lines used in rest-frame 
optical studies. Finally, high-redshift sources are fainter and 
hence more difficult to observe. Probably due to all of these 
possible problems, no reverberation measurements exist for 
AGNs with L > 10 46 ergs s -1 , and several attempts at such 
measurements have so far not been successful (e.g., Welsh 
et al. 2000; Trevese et al. 2004; A. Marconi 2005, private 
communication). 

In view of the many unknowns and the opposing effects en- 
tering the above discussion, and considering the importance 
of the subject, over a decade ago, we began a reverberation- 
mapping program aimed at high-luminosity, high-redshift 
AGNs. In this paper, we describe our program and present 
some initial results. In § 2 we describe the sample, the ob- 
servations, and the data reduction. In § 3 we perform a time- 
series analysis for the AGN light curves and discuss our re- 
sults, with a summary in § 4. Throughout this paper we use 
the standard cosmology with Hq = 70 km s _1 Mpc -1 , Qm = 0.3, 
and fl A = 0.7 (Spergel et al. 2003; Riess et al. 2004). 



2. SAMPLE, OBSERVATIONS, AND DATA REDUCTION 
2.1. Sample Selection 

We selected our sample of 1 1 high-luminosity quasars in 
1994 from the Veron-Cetty & Veron (1993) catalog. The se- 
lection criteria were observationally oriented: a high declina- 
tion (8 > 60°), so as to make the quasars observable from the 
Northern hemisphere for most of the year; observed magni- 
tude V < 18, to permit the acquisition of low -resolution, high 
signal-to-noise (S/N) spectra within a reasonable observing 
time; and, since high-luminosity quasars are known only at 
high redshifts, redshifts in the range 2 < z < 3 .4 that brings the 
prominent broad UV emission lines (C IV and Lyct) into the 
optical range. We chose the 1 1 objects with the brightest ap- 
parent optical magnitudes that satisfy these criteria. The sam- 
ple is listed in Table Q] The choice of bright quasars at high 
redshift naturally favors quasars with high luminosity. The lu- 
minosity range of our sample is 10 464 < ALa(5100) < 10 47 6 
ergs s _1 . This is an order of magnitude higher than other 
AGNs with existing reverberation measurements. 

In FigureQ]we plot the absolute i magnitude versus redshift 
for all objects in the Sloan Digital Sky Survey (SDSS; York et 
al. 2000) Data Release 3 (DR3; Abazajian et al. 2005) Quasar 
Catalog (Schneider et al. 2005), compared to our sample and 
to other objects with reverberation-mapping data. Since the 
SDSS DR3 quasar catalog covers about 10% of the sky, fu- 
ture surveys will not find large numbers of quasars signifi- 
cantly more luminous than those in it, and thus the objects 
in our sample are clearly among the most luminous quasars 
in the Universe. Our sample includes HS 1700+6416 and 
SBS 1425+606 which are the two most-luminous quasars in 
DR3. Our study is probing powerful quasars at redshifts near 
the peak of their comoving number density (near z ~ 2), which 
is the main era of black-hole growth for the most massive 
black holes and galaxies. Figure Q] also shows that there are 
only two objects with M, < -26 that have a previously mea- 
sured BLR size, which illustrates the gap in luminosity be- 
tween our sample and previously studied AGNs. 

The luminosities of the quasars in our sample are likely 
intrinsic, rather than being affected by gravitational-lensing 
magnification. To our knowledge, none of the objects is a 
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FIG. 1. — Absolute i magnitude vs. redshift for all SDSS Data Release 
3 quasars (dots). The six quasars monitored at the HET are marked with 
diamonds, and the five objects from our sample with photometric data only 
are marked with squares. AGNs with reverberation-mapping data are marked 
with triangles. The objects in our sample are clearly among the most lu- 
minous quasars and lie within the primary epoch of black-hole growth for 
luminous quasars. 



lensed quasar. In particular, five of the quasars were im- 
aged with the Hubble Space Telescope (HST) by Maoz et al. 
(1993), who found no evidence for strong lensing. Maoz et 
al. (1993) and other optical lensed-quasar surveys (e.g., Mor- 
gan et al. 2003) have consistently found a lensed fraction of 
~ 1 % among luminous quasars, and thus it is not surprising 
that none among our 1 1 objects is lensed. We also note that 
none of our objects is a BAL quasar, and none has significant 
UV absorption lines. 

Because the Veron-Cetty & Veron (1993) catalog was a lit- 
erature compilation, our sample is not statistically complete, 
nor necessarily representative of the luminous quasar popula- 
tion. For example, six of the quasars in our sample (i.e., about 
55%) are radio loud, somewhat more than the ~ 20% radio- 
loud fraction (i.e., 2 objects) one would expect for a quasar 
population of this mean redshift (2.8) and luminosity (Jiang 
et al. 2006). We note that the incompleteness of AGN rever- 
beration samples applies to most of the studies that have de- 
rived BLR-size-mass-luminosity relations. An exception is 
the sample of 17 PG quasars mapped by Kaspi et al. (2000), 
which were an unbiased selection from the full, flux-limited, 
color-selected Bright Quasar Survey (BQS; Schmidt & Green 
1983; also see Maoz et al. 1994). Of course, the PG sam- 
ple itself suffers from some incompleteness; see Jester et al. 
(2005) and references therein. Although there are no obvious 
biases that are introduced into reverberation measurements by 
these samples (e.g., radio-loud and radio-quiet objects seem 
to follow the same BLR-size-luminosity relation; see the few 
radio-loud objects in the Kaspi et al. 2000 sample and also 
Wu et al. 2004), it should be kept in mind that such biases 
may, nonetheless, exist (see, e.g., Netzer 2003). 

The observations and reductions for the current project fol- 
low the same principles as those of similar projects carried 
out by our group (e.g., Maoz et al. 1994; Netzer et al. 1996; 
Giveon et al. 1999; Kaspi et al. 2000). We outline the main 
features in the following sections. 



2.2. Broad- Band Imaging 

Imaging data have been obtained, starting in 1995, at the 
Wise Observatory (WO; Kaspi et al. 1995) 1 m telescope 
in the Johnson-Cousins B and R bands, using a 1024x1024 
back-illuminated Tektronix CCD. The exposure times are 
250 s in R and 300 s in B. Since the targets have high de- 
clinations, they can be observed from the WO for about 10 
months a year, with observations scheduled about once every 
month (based on the experience with the PG quasars, Kaspi 
et al. 2000, and the expectation that the variability timescales 
of high-luminosity quasars are longer, this sampling rate was 
deemed sufficient). The data have been reduced using IRAF 4 
procedures in the standard way. Broad-band light curves for 
the quasars are produced by comparing their instrumental 
magnitudes to those of constant-flux stars in the field (see, 
e.g., Netzer et al. 1996, for details). The uncertainties on 
the photometric measurements include the fluctuations due to 
photon statistics and the scatter (of 0.03 mag or less) in the 
measurement of the constant-flux stars used. The typical to- 
tal uncertainty is in the range of 0.01-0.05 mag depending 
on object brightness and the observing conditions. The pho- 
tometric light curves provide a check on any continuum vari- 
ability detected in the spectroscopic observations and improve 
the sampling of the continuum light curves. Furthermore, the 
photometric observations began several years before the spec- 
troscopic observations. In the analysis, this measurement of 
the continuum behavior at early times gives an enlarged base- 
line for the cross correlation between the line and continuum 
light curves. 

2.3. Spectrophotometry 

Spectrophotometric monitoring of six of the 1 1 quasars has 
been carried out since 1999 at the 9 m Hobby-Eberly Tele- 
scope (HET; Ramsey et al. 1998; Ramsey et al. in prepa- 
ration). The other five quasars have declinations above the 
HET's high-declination limit (71°41')- The operation of the 
HET in a queue-scheduled mode is a key component of our 
observing program; it enables ideal spreading of the obser- 
vations throughout the observing season for each object and 
ensures that observations will not be lost due to poor weather. 
We have used the Low Resolution Spectrograph (LRS; Hill et 
al. 1998) with a 600 line mm -1 grating, a 10" wide slit, and 
the GG385 order-separation filter. The spectral range covered 
with this configuration is 4300-7300 A with a resolution of 
about 20 A. Typical seeing in the HET data is in the range 
1.5-3". Each night the standard calibration images of bias, 
sky-flat fields, and internal lamp-flats are obtained. Exposures 
of Ne and Cd lamps are used for wavelength calibration (up 
to mid-2002 a Hg-Cd-Zn lamp was used). The objects in our 
sample are accessible to the HET for about six months each 
year. During each such period, we strive to obtain three ob- 
servations separated by ^2.5 months. For each quasar, the 
spectrograph's focal plane is rotated to an appropriate posi- 
tion angle so that a nearby comparison star (columns [9] and 
[10] in Table [TJ is included simultaneously in the slit. This 
permits spectrophotometric calibration of the quasars under 
non-photometric conditions. The wide slit reduces the effects 
of atmospheric dispersion at the nonparallactic position an- 
gles, as well as light losses due to guiding errors and poor 
seeing. 

4 IRAF (Image Reduction and Analysis Facility) is distributed by the Na- 
tional Optical Astronomy Observatories, which are operated by AURA, Inc, 
under cooperative agreement with the National Science Foundation. 
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Observations typically consist of two consecutive expo- 
sures of the quasar/star pair. Exposure times range from 300- 
900 seconds for each exposure depending on object bright- 
ness, aiming for a S/N of > 20 pixel -1 in the continuum. The 
spectroscopic data are reduced using standard IRAF routines. 
The images are bias and flat-field corrected. The extraction 
width is typically 9.4", and wavelength calibration is applied 
to the data after the spectral extraction. This procedure results 
in two spectra from each image: one for the quasar and one for 
the comparison star. The consecutive quasar/star flux ratios 
are then compared to test for systematic errors in the observa- 
tions and to identify and remove cosmic-ray events. The ratio 
is usually reproducible to 0.5^1% at all wavelengths; obser- 
vations with ratios larger than 5% are discarded (this occurred 
in about 7% of the observations). We note a trend of improve- 
ment in the HET observations in image quality and seeing 
over the years due to improvement in its equipment and man- 
ner of operation. We have verified that all our comparison 
stars are non-variable to within by means of differential 
photometry with other stars in each field (see §2.2\ . 

Spectra were calibrated to an absolute flux scale using ob- 
servations of spectrophotometric standard stars observed at 
several epochs. The uncertainty in the calibration is ^20%, 
which is not shown in the error bars in our light curves. The 
error bars reflect only the measurement and differential un- 
certainties, which are of order 1-4%. These are comparable 
to other reverberation mapping studies which have used the 
same observing technique (e.g., Kaspi et al. 1996a; Collier et 
al, 1998; Kaspi et al. 2000). 

For each quasar, we used all available spectra to produce 
an average spectrum and a root-mean-square (rms) spectrum, 
defined as: 



CT < A >= 7^E[/^)-/(A)] 2 , (1) 



where the sum is taken over the N spectra, and /(A) is the 
average spectrum (Peterson et al. 1998a). The average spec- 
tra of the six quasars in our sample are shown in Figure [2] 
We used the average and rms spectra to choose line-free spec- 
tral bands suitable for setting the continuum underlying the 
emission lines, and the wavelength limits for integrating the 
line fluxes. The line and continuum fluxes were measured al- 
gorithmically for all epochs by calculating the mean flux in 
the continuum bands, and summing the flux above a straight 
line in f\ that connects the continuum bands straddling each 
emission line. This process measures the total emission in 
each line, i.e., the flux of the broad component of the line to- 
gether with its narrow component. We note, however, that in 
the high-ionization lines we monitor in luminous quasars, the 
narrow-line contribution is negligible (e.g., Wills et al. 1993). 
In any event, the narrow component is non-variable on the 
timescales probed here, and would merely contribute a con- 
stant flux level to the light curves. 

The uncertainty in the line flux was estimated by propa- 
gating the uncertainty in the definition of the continuum lev- 
els, determined from the standard deviation of the mean in 
the continuum bands (of order 1-3%). To this we added in 
quadrature our estimate for the uncertainty in the calibration 
of the differential spectrophotometry, which is of order 1-4%. 




k . 
5000 6000 7000 

o 

Observed wavelength [A] 



FIG. 2. — Mean spectra of the six quasars monitored with the HET. Each 
spectrum contains roughly a dozen individual observations. The spectral res- 
olution is rj 20 A. The line and continuum bands discussed in the text are 
marked in each panel. 
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TABLE 2 

Integration Limits for Continuum Bands 



Object 


Wavelength Range [A] 


S4 0636+68 


5622-5701 


S5 0836+71 


6456-6556 


SBS 1116+603 


7154-7220 


SBS 1233+594 


6065-6165 


SBS 1425+606 


6051-6141 


HS 1700+6416 


5935-5993 



2.4. Continuum and Emission-Line Light Curves 
As seen in all the quasar spectra in Figure [2] the R band 
(A ~ 6500 A) is dominated by continuum emission, with neg- 
ligible contribution from broad emission lines. In contrast, 
emission lines contribute of order one-half the flux in the 
B band for several of the quasars. The variations in the R- 
band photometry are therefore determined almost completely 
by continuum variations, and this broad-band light curve can 
be merged with the spectroscopic continuum light curve, to 
improve the sampling and duration of the continuum light 
curve. The photometric and spectropho tome trie light curves 
of each object are intercalibrated by comparing all pairs of 
spectrophotometric and photometric observations separated 
in time by less than 20 days. For the data presented here, 
this produced ~ 10 pairs of points per object. A linear least- 
squares fit between the spectrophotometric continuum fluxes 
and the photometric fluxes (the instrumental /^-magnitudes 
were translated to fluxes) was then used to merge the two light 
curves. A clear relation was always apparent in this fitting and 
a linear model was an adequate representation, indicating the 
continuum light curves from the spectroscopic and the photo- 
metric observations are in good agreement. 

For each object, the total number of photometric and spec- 
trophotometric observations is given in Table [T] columns (6) 
& (7), respectively (the sum of the two columns gives the total 
number of points in the continuum light curve). Note that for 
particular emission lines in several objects, some data points 
are missing because of insufficient wavelength coverage or 
low S/N. Light curves for the six quasars monitored at the 
HET are presented in Figure [3] The different bands used for 
the flux measurements in Figure [3] are listed in Table |2] for 
each object (wavelengths are given in the observed frame). 
/?-band light curves for the five objects monitored photomet- 
rically at the WO are presented in Figure H] 

3. TIME SERIES ANALYSIS 

3.1. Continuum and Line Variability 

In Table [3] we list variability measures for all light curves 
presented above. Column (1) lists the object name and col- 
umn (2) the particular Jight curve. In columns (3), (4), and 

(5) we list the mean (/), rms (a), and the mean uncertainty 

(6) of all data points in the light curves, respectively, in the 
appropriate units for each light curve. In column (6) we list 
the xi obtained by fitting a constant to the light curve, and in 
column (7) we list P(x 2 |j/), the probability to get such x 2 if 
there were intrinsically no variability. This test is used to de- 
termine whether the light curve is consistent with a constant 
flux to a significance level of 95%. If the light curve passed 
the test, we list for it in column (8) the intrinsic normalized 
variability measure, an = lOOVc 2 — S 2 // (used by Kaspi et 
al. 2000), and in column (9) we list the fractional variability 



amplitude, F var , and its uncertainty (Rodriguez-Pascual et al. 
1997; Edelson et al. 2002). 

In general, all continuum light curves show variations of 
10-70% measured relative to the minimum flux. The mean in- 
trinsic normalized variability measure (column 8 of Table 
is 5.8% and the median is 4.1%. The points that appear to 
deviate strongly from the mean flux are generally due to un- 
reliable measurements and have relatively large uncertainties. 
A possible exception is one particular deviation in the con- 
tinuum light curve of HS 1700+6416 around JD-2452757 
(2003 April 27) where the object seems to brighten (both in 
the spectroscopic and in the photometric measurements) by 
about 10% for ~ 15 days. This object was monitored also by 
Reimers et al. (2005), who found that on 2003 May 14-18 
HS 1700+6416 was 20% brighter in the R band than in March 
1995 and July 1998. This agrees with our measurements at 
the same epochs and suggests that the above deviation in the 
light curve is real. 

It is instructive to compare the variability characteristics 
of our sample to those of the low-redshift PG quasar sample 
(Kaspi et al. 2000). However, one needs to consider the dif- 
ferences in the rest-frame monitoring periods of the two sam- 
ples. While the continuum light curves of the high-luminosity 
quasars span over 10 yr in the observed frame, this period 
corresponds to only about 2-3 yr in their rest frames. We 
have therefore analyzed the first 3 yr period in the PG quasar 
light curves (which span over 7.5 yr). The resulting variabil- 
ity measured for the PG quasar sample, 25-125%, is about 
double the 10-70% variability range of the high-luminosity 
quasars. Also, the mean and rms of for the PG quasars (in 
the first 3 years) is <aAr>=(10.4±4.5)% about twice that of 
the high-luminosity quasars which have <(7;v>=(5.8±3.9)%. 
The lower rest-frame variability we measure in the continuum 
for the current sample is probably a manifestation of the gen- 
eral trend that high-luminosity AGNs have longer variability 
timescales (e.g., Vanden Berk et al. 2004), perhaps as a result 
of their higher black hole masses. 

Proceeding to the emission-line light curves, the variabil- 
ity measures listed in Table [3] show that none of the four Lya 
light curves shows significant variability. In contrast, the two 
C III] A 1909 s light curves and all six C IV light curves show 
significant variability, with F var significant at the >3er level. 
We note that the variability measures of these emission-line 
light curves are comparable to, or even greater than, those 
of their corresponding continuum light curves. There are few 
previous AGN UV data sets with which to compare these pos- 
sible trends. The only quasar with UV variability data of 
similar quality is 3C273. Interestingly, Ulrich et al. (1993) 
noted the non-variability of Lya in this object, at a level of 
< 5%, over a period of 15 years, despite factor-of-two vari- 
ations in the continuum during the same period (there are no 
data for the C IV and C III] lines during that time). In contrast, 
Seyfert galaxies seem to show comparable variation ampli- 
tudes in the UV continuum and in the Lya, C IV, and C III] 
lines, though the line-variation amplitudes tend to be some- 
what smaller than the continuum variations. This has been 
seen in NGC 5548 (Clavel et al. 1991; Korista et al. 1995), 
NGC 7469 (Wanders et al. 1997), NGC 3783 (Reichert et al. 
1994), Fairall 9 (Rodriguez-Pascual et al. 1997) and 3C390.3 
(O'Brien et al. 1998). 

The indication that, in quasars, the C IV, C III], and Lya 

5 The C III] complex we measure includes up to ~50% flux contribution 
from Si III]A1892 and Fe III multiplets. 
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FIG. 3. — Light curves for the six quasars which have been monitored at the HET. Squares are spectrophotometric data from the HET. Triangles are photometric 
data from WO. Time is given in Julian Day (bottom) and UT date (top). Continuum flux densities, fx, are given in units of 1CT 16 ergs citT 2 s _1 A -1 and emission- 
line fluxes are given in units of 1CT 14 ergs cm -2 s -1 . 
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FIG. 4. — R-band light curves for the five quasars which have been mon- 
itored photometrically at the WO. Magnitudes are given in arbitrary units. 
Time is given in Julian Day (bottom) and UT date (top). 



lines show different variability trends may be related to the 
fact that each line is driven by a different ionizing continuum. 
This could also explain the different variability amplitudes, 
e.g., if the Lya line is driven by a weakly varying UV con- 
tinuum, while the C IV is controlled by a more strongly vary- 
ing extreme-UV to X-ray continuum. Indeed, the continuum 
variability amplitudes of AGN are observed to grow with de- 
creasing wavelength, from optical through the UV and to the 
X-rays, and the variations are sometimes poorly correlated 
with each other in different bands (e.g., Maoz et al. 2002; 
Shemmer et al. 2003). The differences in response behaviors 
among quasars and Seyferts would then imply that the degree 
to which continuum variations grow with decreasing wave- 
length depends on luminosity. However, the current data are 
still sparse, and the above trends require confirmation. 

Comparison of the variability in the emission-line light 
curves between the PG sample and the present high- 
luminosity sample is complicated by the fact that, in the PG 
sample, the Balmer lines were monitored, while for the high- 
luminosity quasars we monitor UV lines. Furthermore, the 
high-luminosity quasar lines were monitored for only half the 
period over which the continuum was monitored, i.e., approx- 
imately 5 yr, which is about 1.2-1.5 yr in the quasar rest 
frames. Analyzing only the first 1 .5 yr of the Balmer-line light 
curves of the PG quasars, we find that about one-third of the 
lines did not vary; most of the non-varying light curves are for 
the H7 line which has lower EW than the Ha and H/3 lines. 
If, for the PG-quasar emission lines, we average the variabil- 
ity measures of only Ha and H/3, we find <<tn>=(7 .5±6.2)% 
(mean and rms), compared to the average variability measure 
of the C III] and C IV lines of the high-luminosity quasars in 
the current sample, <otv>=(12.8±9.4)%. It therefore appears 
that the UV lines may be more responsive than the Balmer 
lines to continuum variations, which could compensate for the 
lower continuum variability amplitudes at high luminosities, 



and contribute to the success of our experiment. Indeed, the 
high ratio of emission-line to continuum-variation amplitudes 
in the high-luminosity sample is in contrast to the common 
behavior in the PG quasars which show less variability in the 
Balmer lines than in the continuum. [Some low-luminosity 
AGNs, e.g., NGC 405 1 (Shemmer et al. 2003), do show larger 
variability amplitudes in the Balmer lines than in the optical 
continuum.] However, due to the limited rest-frame period 
considered here, we consider this result tentative. 

We do not see any strong FWHM or line-profile variations 
in our sample of high-luminosity quasars. However, at the low 
resolution of our spectra (about 1200 kms -1 at 5000 A), our 
ability to detect such variations is limited to dramatic changes 
only. 

3.2. Line -Continuum Cross Correlations 

The ultimate objective of our program is to detect and mea- 
sure a time delay between the continuum and the line-flux 
variations in high-luminosity AGNs. The significant contin- 
uum and line variations that we have observed during a decade 
demonstrate that, at least in principle, such a measurement 
may be feasible. However, an actual estimate of a delay re- 
quires: (1) high variability amplitude in both the line and 
the continuum light curves; and (2) light curves with at least 
a few, reasonably sampled, large-amplitude "events", i.e., 
changes of sign in slope, that permit unambiguous matching 
between continuum and line variations. Examining the light 
curves of the six quasars with emission-line data at the cur- 
rent stage of our project, all but one currently suffer from ei- 
ther low variability amplitude in the emission-line light curves 
or monotonically increasing or decreasing continuum light 
curves. In four objects the C IV light curve tracks the contin- 
uum's monotonic variation (SBS 1 1 16+603, SBS 1233+594, 
SBS 1425+606, HS 1700+6416). We note that the monotonic 
variation trend for these objects is unlikely to continue indef- 
initely. We expect that at some point, there will be a "break" 
in each of these cases, permitting the eventual measurement 
of a lag. However, the one current exception to this state of 
affairs is S5 0836+71; although the data for this quasar are 
still not ideal in terms of the criteria above, they do allow a 
preliminary measurement of the emission-line to continuum 
lag. The continuum light curve for this object displays a gen- 
eral rise until June 2002, followed by a sharp drop in flux. 
This light curve has largest variation among all our monitored 
quasars, and the light curves of both C IV and C III] display 
the largest variations among all the emission line light curves. 
Both emission-line light curves seem to follow the general 
trend of the continuum light curves, although the C IV light 
curve seems to have a much larger time lag than the C III] 
light curve. 

To quantify the time lag, we use two methods for cross- 
correlating the line and continuum light curves. The first 
method is the interpolated cross-correlation function (ICCF), 
as implemented by White & Peterson (1994; see also the 
review by Gaskell 1994). The second method is the z- 
transformed discrete correlation function (ZDCF) of Alexan- 
der (1997). The two methods yield similar results for the cur- 
rent data, and we will use only the ICCF results in the follow- 
ing analyses. The results of the cross-correlation analysis are 
presented in Figure [5] 

The CCFs for S5 0836+71 show enough structure to permit 
estimation of a time lag. The clearest case is the correlation 
between the C IV line and the continuum. To quantify the 
possible time lag and its uncertainty, we use the the model- 
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TABLE 3 
Variability Measures 
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NOTE. — No data in columns (8) & (9) mean that no significant variability was detected in this light curve. In particular none of the Lya 
light curves shows significant variability 

d The units are 10" ergs cm - s" A" for the continuum light curves, 10" ergs cm"" s" for the emission-ine light curves, and instrumental 

magnitude for the /?-band light curves. 

k Probabilities smaller than 10" 9 are listed as zero. 

independent FR/RSS Monte Carlo method of Peterson et al. 
(1998b; 2005). In this method, each Monte Carlo simulation 
is composed of two parts. The first is a "random subset selec- 
tion" (RSS) procedure which consists of randomly drawing, 
with replacement, from a light curve of TV points a new sample 
of N points, while preserving the temporal order. The second 
part is "flux randomization" (FR) in which the observed fluxes 
are altered by random Gaussian deviates scaled to the uncer- 
tainty ascribed to each point. The two resampled and altered 
time series are cross-correlated using the ICCF method, and 
the centroid of the CCF is computed. We used ^10000 Monte 
Carlo realizations to build up a cross-correlation centroid dis- 
tribution (CCCD; e.g., Maoz & Netzer 1989). The mean of 
the distribution is taken to be the time lag and the uncertainty 
is determined as the range which contains 68% of the Monte 
Carlo realizations in the CCCD, and thus would correspond 
to 1(7 uncertainties for a normal distribution. 

We find the time lag between the C IV line and the con- 
tinuum of S5 0836+71 to be 595+ff days, or 188+^ days in 
the quasar rest frame. However, 17% of the CCCD simu- 
lations failed to produce significant centroid measurements. 
This may indicate that the quality of the data render this anal- 
ysis premature. The cross correlation of the C III] light curve 
with the continuum yields a peak around zero time lag. A 
formal CCCD gives -152l[|| days, which is -48+^ days in 
the rest frame, but the CCCD function is not simple, with sev- 
eral peaks within this uncertainty range. We note that the only 
AGNs in which (rather uncertain) C III] time lags have been 
estimated to date are NGC 5548 (Peterson et al. 2004, and 
references therein) and NGC 4151 (Metzroth et al. 2006). 

Until recently, only four AGNs had measured C IV rever- 
beration time lags: NGC 3783, NGC 5548, NGC 7469, and 
3C 390.3 (see Peterson et al. 2004, for a summary). Peterson 
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FIG. 5. — Cross-correlation functions, ICCF (solid curve) and ZDCF (cir- 
cles with error bars), between the continuum and the emission-line (C IV - 
top, C III] - bottom) light curves of S5 0836+71 from Fig.[5] 
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FIG. 6. — Size-luminosity relationship based on the C IV A 1549 emission 
line and the UV continuum. Data (circles) are reproduced from Peterson et al. 
(2005), to which we add S5 0836+71 from the current study. The new data 
point deviates somewhat from the best-fit relation found by Peterson et al. 
(dotted line). Our own linear fits to the data are marked with a solid line for 
the FITEXY method and a dashed line for the BCES method. The power-law 
slopes of the three relations are indicated. 



et al. (2005) have now measured the C IV time lag also for the 
low-luminosity Seyfert galaxy NGC4395, which is four or- 
ders of magnitude lower in luminosity than those four AGNs, 
permitting a first estimate of the BLR radius-luminosity re- 
lationship for the C IV-emitting region. Our preliminary de- 
termination of the C IV time lag for S5 0836+71 allows us 
to extend this relation to seven orders of magnitude in lumi- 
nosity. In Figure [6] we show the data as presented by Pe- 
terson et al. (2005), as well as their adopted best-fit rela- 
tion (dotted line), to which we have added the data point for 
S5 0836+71. The rest-frame UV luminosity of S5 0836+71 
has been computed from the mean flux during our observa- 
tions. We have corrected the luminosity for a Galactic ex- 
tinction of A v = 0.101 from Schlegel et al. (1998) using the 
extinction curve of Cardelli et al. (1989). We find the lumi- 
nosity to be AL A (1350A)= (1.12±0.16) x 10 47 ergs s -1 . Our 
new data point deviates from the extrapolation of the fit by 
Peterson et al. (2005), which would predict a rest-frame time 
delay for S5 0836+71 of 372-1865 days. We note that, due 
to the current length of our program, we are sensitive to time 
lags of up to about 1000 days in the observer frame. 

Most previous reverberation-mapping experiments have 
measured the Balmer emission lines, in particular H/3. For 
the four AGNs above, H/3 reverberation measurements exist, 
in addition to the C IV measurements (Peterson et al. 2004). 
The scaling between the BLR sizes of the two lines is im- 
portant since it is the means by which the BLR size from 
C IV reverberation mapping of high-luminosity quasars can 
be compared to the BLR size from H/3 reverberation map- 
ping of Seyferts and low-luminosity quasars. For NGC 3783, 
NGC 5548, and NGC 7469 the H/3 delays are factors of 2.55, 
1.98, and 1.8, respectively (which average to ~ 2) larger than 
the C IV delays, but for 3C 390.3 the C IV delay is larger than 
the H/3 delay by a factor of 1 .5. Thus, the relation between the 
C IV and H/3 BLR sizes is still unclear, due to the small num- 
ber of objects with reverberation measurements in both lines, 
and the small spread in luminosity among these four AGNs. 
BLR stratification suggests that the C IV emission should arise 
interior to the H/3 emission. This is confirmed by the three 



Seyfert Is but not by the measurements of 3C 390.3. Also, 
the very small spread of these four AGNs in luminosity does 
not allow a meaningful scaling of the H/3 BLR size to the C IV 
BLR size; broadening the luminosity range is essential. 

We have performed a linear regression for the C IV BLR- 
size-luminosity relation using all current data points in Fig- 
ure [6] We have used the two methods which are described in 
detail in Kaspi et al. (2005): (1) the FITEXY from Press et al. 
(1992, p. 660) which is based on an iterative process to min- 
imize x 2 an d m which we adopted the Tremaine et al. (2002) 
procedure to account for the intrinsic scatter in the data; (2) 
the bivariate correlated errors and intrinsic scatter (BCES) re- 
gression method of Akritas & Bershady (1996). Both meth- 
ods take into account the uncertainty in both coordinates and 
the intrinsic scatter around a straight line. The fit for the data 
we find from the FITEXY method (solid line in Figure [6) is 

^^=(0.17 ± 0.04)f^H^V 52±a ° 4 , (2) 
lOltdays V 7 V 10 ergss" 1 J ' 

and for the BCES method (dashed line in Figure|6]) it is 

0.55±0.04 



= (0.24 ± 0.06)f AL 4 f 350A 1 ) ) 
1 10 43 ergss" 1 J 



(3) 

lOltdays " " V 10 43 ergss"' ' 

The Pearson correlation coefficient is 0.97 with significance 
level of 7 x 10~ 5 and the Spearman rank-order correlation is 
0.83 with a significance level of 0.01. The intrinsic scatter 
we find using the Tremaine et al. (2002) procedure is 54%. 
The slope we find for the relation between the C IV BLR-size 
and UV luminosity is very close to the one found by Kaspi 
et al. (2005) between H/3 BLR-size and the UV luminosity 
(-0.55). 

The mean FWHM of the C IV line measured from the mean 
spectrum of S5 0836+71 is about 9700 km s" 1 . Using Eq. 5 of 
Kaspi et al. (2000) and the time lag of 188 days, we estimate 
the central mass of S5 0836+71 at about 2.6 x 10 9 M o . This 
is the highest mass directly measured for the black hole in an 
AGN using reverberation mapping. 3C273 (=PG 1226+023), 
the quasar with the highest directly measured mass so far, has 
a mass of 8.9 x 10 8 M Q , AL A (1350 A)=2.0x 10 46 ergs s" 1 , and 
AL A (5 100 A)=9. 1 x 10 45 ergs s" 1 . Thus, S5 0836+71 has a fac- 
tor three higher mass and a factor — 6 higher luminosity than 
3C273. 

The central mass of S5 0836+71 corresponds to an Ed- 
dington luminosity of Lgdd = 3.9 x 10 47 ergs s" 1 (for solar- 
abundance gas L Edd w 1.5 x 10 38 M/M Q , e.g., Shaviv 1998). 
Using AL A (5100 A)=6.4x 10 46 ergs s" 1 for S5 0836+71 and a 
bolometric correction of 5.6 from Marconi et al. (2004, Eq. 
21) we find a bolometric luminosity, Lbol = 3.6 x 10 47 ergs 
s . Thus, this object is accreting close to its Eddington lu- 
minosity, with Lboi/^Edd ~ 0.9. We note that using Eq. 7 
from Vestergaard & Peterson (2006) to estimate the black- 
hole mass from the UV luminosity and the C IV FWHM, one 
obtains Mbh = 1.8 x 10 10 M Q , which is 7 times higher than 
our tentative measurement from reverberation mapping. 

The high Eddington ratio that we obtain is similar to that 
found for many other quasars (e.g., Shemmer et al. 2004; 
Baskin & Laor 2005). However, we note that there are several 
caveats which could influence our result. The first is the ten- 
tative nature of our measurement of the time lag for the C IV 
line. One should keep in mind the possibility that the corre- 
lation and delay that we have measured in S5 0836+71 might 
be a chance coincidence, since it is the single object in which 
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we have found a correlation, out of six objects where we have 
searched for one. Another concern is our measured luminos- 
ity of the object. Since S5 0836+71 is the radio-loudest object 
in our spectrophotometric sample it could be that a jet-linked, 
beamed optical continuum is influencing the variability and 
the continuum-flux measurement. While a beamed compo- 
nent might influence the variability on small timescales (the 
continuum light curve does show variability on timescales of 
months) the C IV line light curve is smooth and its response is 
to the long-timescale variations seen in the continuum light 
curve. We have also checked the C IV equivalent widths 
(EWs) of the six quasars monitored at the HET, of which 
three are radio loud and three are radio quiet. For the radio- 
loud quasars, the observed EWs are 33, 68, and 101 A while 
for the radio-quiet ones the EWs are 57, 68, and 90 A. Thus, 
both classes have comparable C IV EWs, and there is no clear 
difference that might indicate that there is an additional con- 
tribution to the UV continuum luminosity of the radio-loud 
quasars. 

4. SUMMARY 

We have presented first results from our long-term program 
to monitor high-luminosity, high-redshift quasars, with the 
objective of eventually measuring the sizes of their broad-line 
regions and the masses of their central black holes. All 1 1 
quasars in our photometrically monitored sample show con- 
tinuum variability of 10-70% over a rest- frame period of 2- 
3 years. Compared to previous studies of lower luminosity 
AGNs, the continuum variability amplitude of higher lumi- 
nosity quasars over the same rest-frame interval appears to 
be smaller in the current high-luminosity sample. Six of the 
1 1 quasars have been monitored spectroscopically over the 
past ^5 years (1.2-1.5 yr in the quasar rest frames) using the 
HET. None of the quasars show Lya variations, at a level of 
< 7%, and this non-variability may be a generic feature of 
high-luminosity AGN. However, we do detect variable C III] 
and C IV broad line fluxes, whenever these lines are in our 
observed spectral range. In several cases, these lines track 
the continuum variations in the same quasar. The variations 
in the broad C III] and C IV lines are higher than Balmer-line 
variations in low-luminosity quasars over the same rest-frame 
interval, but this result is tentative, in view of the still-limited 
spectroscopic monitoring period. 

In one object, S5 0836+71, we tentatively measure the C IV 
rest-frame time lag behind the continuum emission to be 
1 88+37 days. This is lower than expected from an extrap- 
olation of the C IV BLR-size-luminosity relation measured 
at lower luminosities. Our measurement permits tracing this 
relation over seven orders of magnitude in AGN luminosity. 
From this time lag we estimate a black hole mass of 2.6 x 10 9 
Mq in S5 0836+71, the largest ever measured by this tech- 
nique. 

Our results to date demonstrate that reverberation mapping 



of high-redshift, high-luminosity, quasars may be feasible. 
Based on our past experience at quasar reverberation (Maoz 
et al. 1994; Kaspi et al. 1996b, 2000), we hope that with 
ss 5 more years of similar data we will be able to establish 
a reliable BLR size for S5 0836+71, and perhaps for several 
other objects in our sample. The unpredictability of quasar 
light curves, together with our lack of knowledge about the 
response amplitudes and timescales of the UV lines, makes 
it difficult to estimate our chances of success. If we do suc- 
ceed, the combined reverberation results from this program 
and those from a program being carried out on very low- 
luminosity AGNs will cover the entire AGN luminosity range 
and much of cosmic time, including the peak era of black- 
hole growth for the most massive black holes and galaxies. 
This will enable a direct measurement of the black-hole mass 
in these AGNs and a more reliable indirect mass measurement 
for all AGNs. 
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